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ABSTRACT 

A spatial clustering signal has been established in Spitzer/IRAC measurements of the unresolved cosmic near- 
infrared background (CIB) out to large angular scales, ~1°. This CIB signal, while significantly exceeding the 
contribution from the remaining known galaxies, was further found to be coherent at a highly statistically significant 
level with the unresolved soft cosmic X-ray background (CXB). This measurement probes the unresolved CXB to 
very faint source levels using deep near-IR source subtraction. We study contributions from extragalactic populations 
at low to intermediate redshifts to the measured positive cross-power signal of the CIB fluctuations with the CXB. 
We model the X-ray emission from active galactic nuclei (AGNs), normal galaxies, and hot gas residing in 
virialized structures, calculating their CXB contribution including their spatial coherence with all infrared emitting 
counterparts. We use a halo model framework to calculate the auto and cross-power spectra of the unresolved 
fluctuations based on the latest constraints of the halo occupation distribution and the biasing of AGNs, galaxies, 
and diffuse emission. At small angular scales (<F), the 4.5 /im versus 0.5-2 keV coherence can be explained by 
shot noise from galaxies and AGNs. However, at large angular scales (~10 / ), we find that the net contribution 
from the modeled populations is only able to account for ~3% of the measured CIB x CXB cross-power. The 
discrepancy suggests that the CIB x CXB signal originates from the same unknown source population producing 
the CIB clustering signal out to ~1°. 

Key words: cosmic background radiation - dark ages, reionization, first stars - diffuse radiation - 
large-scale structure of universe 
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1. INTRODUCTION 

Background radiation is produced by a variety of sources 

that dominate different parts of the electromagnetic spectrum 

and includes contributions from sources inaccessible to indi- 
vidual telescopic studies. The near-IR cosmic infrared back- 
ground (CIB) probes the history of starlight and associated emis- 
sion falling into the 1-5 /im range (see review by Kashlinsky 
2005), whereas the cosmic X-ray background (CXB) probes 
emissions from accretion-powered sources and hot ionized gas. 
Spatial correlations between the two can arise from sources 
emitting at both IR and X-ray wavelengths or from sepa- 
rate IR and X-ray sources that share the same large-scale 

structures. Angular fluctuations in the CIB have been re- 
vealed after carefully masking resolved sources in deep near- 
infrared (near-IR) exposures and Fourier analyzing the pixels 
remaining in the maps. This technique has led to measure- 
ments using Spitzer/IRAC, Hubble Space Telescope / N I CMOS , 
and AXARZ/InfraRed Camera (Kashlinsky et al. 2005, 2007a; 
Thompson et al. 2007; Matsumoto et al. 2011; Kashlinsky et al. 
2012; Cooray et al. 2012b). The amplitude of the CIB fluctu- 
ations in the Spitzer and AKARI data, now extending to ~1°, 
is well above the expected contribution of local foregrounds, 

such as the Zodiacal Light and Galactic cirrus, and implies 
an unresolved extragalactic CIB component which significantly 
exceeds the power from the remaining known galaxy popula- 
tions (Helgason et al. 2012). The minimal CIB flux implied by 
the new sources responsible for these fluctuations is of order 


0.5-2 nW in 2 sr -1 (Kashlinsky et al. 2007a) and is well below 
the claimed direct CIB flux measurements from DIRBE (Dwek 
& Arendt 1998), Infrared Telescope in Space (Matsumoto et al. 
2005), and AKARI (Tsumura et al. 2013), being consistent 
with limits placed by y-ray attenuation from very high energy 
sources (Meyer et al. 2012; Ackermann et al. 2012; H.E.S.S. 
Collaboration 2013). The power spectrum of the Spitzer fluc- 
tuations is highly isotropic on the sky and is consistent with 
that of high- 7 ACDM clustering out to ~1° (Kashlinsky et al. 
2012); it appears, within the statistical uncertainties, with the 
same amplitude and spatial shape at both deeper (lower shot- 
noise; Kashlinsky et al. 2007b) and shallower (larger shot-noise; 
Cooray et al. 2012b) IRAC maps. 

The measured spatial power spectrum of the source- 
subtracted CIB fluctuations rises at angular scales >20". The 
amplitude and shape of this rise is a direct measurement of 
the clustering properties of the underlying source populations 
and thus a primary key of understanding their nature. The in- 
terpretation of these power spectra has been a matter of de- 
bate (Salvaterra et al. 2006; Kashlinsky et al. 2007b; Fernandez 
et al. 2010; Cooray et al. 2012a; Yue et al. 2013a). It is now 
firmly established that the extragalactic signal is inconsistent 
with the emission of presently resolved galaxies (z < 6) with 
extrapolation to low luminosities (Helgason et al. 2012). On 
the other hand, it was proposed that CIB fluctuations from 
the first light era could be measurable thus making the sig- 
nal a critical tool for studying the high-7 universe (Kashlinsky 
et al. 2004; Cooray et al. 2004). The CIB fluctuations show 
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no spatial correlations with deep Hubble / Advanced Cam- 
era for Surveys (ACS) maps <0.9 p m implying that the 
sources are extremely faint and/or that their Lyman-break 
is redshifted well into the near-IR, (1 + z)0. 1 p m > 0.9 /im 
(Kashlinsky et al. 2007c). Current measurements show strong 
clustering on large scales coupled with the low shot noise lev- 
els on small scales (Kashlinsky et al. 2007b). This is consistent 
with a high redshift origin, although an alternative scenario has 
recently been proposed to explain some of these measurements 
as arising from the intrahalo light from stars stripped of their 
paternal halos at intermediate redshifts of z ~ (1—4) (Cooray 
et al. 2012b; see Section 4.3 for full discussion). 

The large-scale clustering of the source-subtracted CIB fluc- 
tuations does not provide direct information on whether the 
underlying sources are powered by stellar nucleosynthesis or 
accretion onto compact objects. Cappelluti et al. (2013, here- 
after Cl 3) provided observational evidence for a substantial 
population of accreting sources among the CIB sources, rais- 
ing the intriguing possibility of extensive black hole activ- 
ity in the early Universe. C13 used deep source-subtracted 
Spitzer/IRAC and Chandra maps of a common region to re- 
veal a highly significant cross-power (3.8a, 5.6a) between the 
unresolved CIB at 3.6, 4.5 pm and the soft 0.5-2 keV CXB, 
arguing for a high-z origin of the sources. An interesting and 
specific model for the discovered CIB x CXB signal explained 
it in terms of direct collapse black holes at z > 12-15 (Yue 
et al. 2013b) by a mechanism which is capable of reproduc- 
ing both the unexplained CIB and cross CIB x CXB fluctuation 
signals without violating constraints imposed by the total mea- 
sured soft CXB. However, before models of such hypothetical 
high-z sources can be favored as leading explanations for the 
measured signal, a more quantitative analysis of known, and 
proposed, source classes at z < 6 is needed. There exists a vari- 
ety of known mechanisms of X-ray production which can spa- 
tially correlate with optical/IR emitting counterparts in complex 
ways. 

The deepest Chandra surveys have been able to resolve 
~80%-90% of the [0.5-8] keV CXB into individual point 
sources, the majority of which is made up of active galactic 
nuclei (AGNs; Hickox & Markevitch 2006; Lehmer et al. 
2012). However, at the faintest fluxes, the abundance of sources 
identified as normal galaxies rapidly approaches that of AGNs 
and is likely to dominate at fainter levels. The X-ray emission 
within galaxies comes predominantly from X-ray binaries 
(XRBs), a compact object accreting from a companion star, 
which have been found to scale well with galaxy properties, 
such as star formation rate (SFR) and stellar mass (e.g., 
Ranalli et al. 2003; Persic & Rephaeli 2007; Lehmer et al. 
2010). These sources have recently been detected out to deeper 
levels and higher redshifts by stacking analyses (Cowie et al. 
2012; Basu-Zych et al. 2013). Cappelluti et al. (2012b) studied 
the unresolved CXB fluctuations (>2 x 10~ 16 erg s -1 cm -2 ) 
remaining in deep Chandra exposures and determined that 
the bulk is produced by gas residing in galaxy groups and 
clusters (~50%), with the rest being contributed by AGNs 
and galaxies. At these levels, any contribution from high-z 
miniquasars would be overwhelmed by these low-z components 
(<5%), although the systematic uncertainty in the mean level 
of the CXB increases such constraints by a large factor. The 
signal measured in C13 is revealed only after eliminating 
undetected X-ray sources down to unprecedented flux levels, 
<JC5 x 10” 17 erg s ' 1 ctrT 2 , using the deep source-subtracted 
Spitzer maps. 


In this paper, we explore the contribution of the intermediate 
Z sources to the measured level of the CIB x CXB coherence 
(Cl 3) by modeling components from different populations: 
galaxies, AGNs, and diffuse emission in both IR and X-rays. 
We also discuss the intrahalo model of Cooray et al. (2012b) 
versus the full set of CIB constraints and its measured coherence 
with the unresolved CXB. This exhausts the set of known and 
proposed populations out to z < 4-6. We present a formalism 
for reconstructing the cross-power spectrum of the fluctuations 
produced by each source class using the latest observational 
evidence for their clustering and abundance. 

The paper is organized as follows. Following Section 2, 
which defines the basic parameters, Section 3 discusses the 
CIB xCXB measurement of C13 in more detail. Section 4 
addresses contributions to the CIB and presents population 
modeling of optical/IR emitting sources. In Section 5, we do 
the same for CXB contributions with X-ray population models. 
In Section 6, we develop the formalism for reconstructing 
the CIB xCXB fluctuation signal and present our results in 
Section 7, followed by discussion in Section 8. 

2. DEFINITIONS AND PARAMETERS 

Throughout, all quantities referred to as “X-ray” or denoted 
by “X” correspond to the emission in the soft X-ray band 
0.5-2 keV, unless noted otherwise. The suffix IR refers to 
the near-IR wavelengths 3.6 pm and 4.5 p m , i.e., the effective 
wavelengths of Spitzer/IRAC bandpasses 1 and 2, whereas 
FIR refers to the total infrared quantities, integrated over 
10-1000 pm. We assume standard ACDM cosmology using the 
cosmological parameters from the Planck experiment (Planck 
Collaboration et al. 2013). All magnitudes are in the AB system. 

Background fluctuations are characterized by the spatial 
power spectrum as a function of the spatial frequency q (or 
spatial scale 2tz / q), defined as P(q ) = (|A(</)| 2 ), where A(q) is 
the two-dimensional Fourier transform of the source-subtracted 
CIB. The mean square fluctuation in CIB flux on angular 
scale 6 = 2 n/q is defined as (&F 2 ) = q 2 P(q)/(/bx). The 
cross-power describing the correlations between fluctuations 
at different wavelengths ( m,n ) is P m „(q ) = (A m {q)A*(q)) = 
1Z m {q)1Z n {q) + l m (q)T n (q) with 1Z , I standing for the real, 
imaginary parts of the Fourier transform, A{q). The cross- 
power spectrum is a real quantity which can assume positive 
or negative values. Coherence is then defined in its usual way as 
in Kashlinsky et al. (2012), C{q ) = ([ P mn {qj] 2 / P m (q)P n {q ))• In 
the absence of common (coherent) populations at wavelengths 
m and n, the cross-power, measured from a map of /V p j x pixels, 
will oscillate around zero with a random statistically uncertainty 
of order [P m P n ]/ pix- Table 1 summarizes the definitions of 
the quantities used throughout the paper. 

The source-subtracted CIB fluctuations have two compo- 
nents: (1) large scales >(20"-30") are dominated by the clus- 
tering of an unresolved population of unknown sources, while 
(2) small scales are dominated by the shot(white) noise, 7 which 
arises mainly from unresolved galaxy populations. The two 
components may arise from different populations of unresolved 
sources. 

3. THE MEASURED CIB xCXB COHERENCE 

The spatial coherence measured between the source- 
subtracted CIB fluctuations and the unresolved CXB used data 

7 This refers to the net white noise component, which includes the counting 
noise and, in some models, also the one-halo component. 
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Name 

Symbol 

Expression 

Units 2 

Infrared flux b 

F 

vl v 

nW m -2 sr -1 

X-ray flux c 

S 

Ie ? N(E)dE 

—1 -2 —1 
erg s cm sr 

Angular scale 

e 

2 /q 

arcsec (") 

Angular wavenumber 

q 

kd c 

rad -1 

Sky brightness 

F(d) 

(F)+8g 

nW m -2 sr -1 

Two-dimensional Fourier transform 

\ 

f 8ge~ iq ' e d6 

nW m” 2 rad -1 

Angular power spectrum 

P(q) 

<AA*> 

nW 2 m” 4 sr” 1 

Cross-power spectrum 

P ) 


erg s _1 cm -2 nW m -2 sr -1 

Fluctuations 

<**»> 

q 2 P(q)/27Z 

nW 2 m” 4 sr” 2 

Broad band averaged power 

(P) 

fql P(.q)qdq / /« 2 qdq 

erg s -1 cm -2 nW m~ 2 sr -1 

Coherence 

C(q) 

P 2 n (q)/(Pm(.q)Pn(q )) 



Notes. 

2 These are the units we use throughout the paper. 
b Refers to Spitzer / IRAC 3.6 /xm and 4.5 /am flux. 
c Refers to soft X-ray 0.5-2 keV band flux. 


from the deep Chandra ACIS-I AEGIS-XD survey and the 
Spitzer Extended Deep Survey (SEDS) in the EGS field, where 
the two data sets overlap in a ~8' x 45' region of the sky (for 
details, see Cappelluti et al. 2013). The measured cross-power 
between the source- subtracted IRAC maps at (3. 6,4. 5) /xm and 
Chandra [0.5-2] keV maps was detected at angular scales of 
10"-1000", with an overall significance of ~(3.8a, 5.6a), re- 
spectively. At the same time, no significant signal was de- 
tected between the IRAC source-subtracted maps and the harder 
Chandra bands. The measured coherence signal has been de- 
tected after jointly masking resolved sources down to m AH ~ 25 
and 7 x 1 O' 1 7 erg s -1 cm~ 2 in IR and X-rays, respectively. The 
signal is characterized by the cross-power spectrum, Pir.xO?) 
shown in Figure 10, and exhibits a broadband averaged cross- 
power in the 10"- 1000" angular range (see definitions in Table 1) 

(F > 3.6 |[i m,0.5-2kev) = 6.4 ± 1.7, 

(F , 4.5 |[im ,o.5-2kev) = 7.3 ± 1.3, 

in units of 10~ 20 erg s -1 cm -2 nW m -2 sr” 1 . We refer to the 
opposite ends of the measured angular interval [10", 1000"] as 
small and large angular scales, respectively. 

In this paper, we examine the contribution to this signal from 
extragalactic populations at z < 6 with the putative high-z 
populations deferred to a forthcoming paper. We decompose 
the total power spectrum of the fluctuations into the sum of 
power from sources of types known to emit both in X-rays and 
optical/IR 

Gotal = Fgaiaxies + F*AGN ^diffuse (1) 

with each of these components contributing both in terms of 
their large-scale clustering on the sky as well as shot noise 
dominating small-scale power. We refer to sources of X-ray 
emission arising collectively from stars, stellar remnants, and 
gas within galaxies as “normal galaxies” whereas the term 
“AGN” is used in its broadest sense referring to any black hole 
activity in the centers of galaxies regardless of subclasses. We 
also consider nonpoint sources, such as hot diffuse gas and 
dispersed light around galaxies, which we collectively refer to 
as “diffuse” components. 

The measured coherence can be interpreted as the product of 
the fraction of the emission due to common populations, i.e., 


C = £ 2 £ 2 , where is the fraction contributed by the common 
populations in band m. At large angular scales (>20"), the level 
of the measured coherence between 4.5 /xm and [0.5-2] keV 
is measured to be C(q) ~ 0.03-0.05, implying that at least 
15%-25% of the large scale power of the CIB fluctuations 
is correlated with the spatial power spectrum of the X-ray 
fluctuations, i.e., -~/C > 15%-25%. This implies that the true 
nature of the signal lies somewhere in between two limiting 
scenarios: 

1. 100% of the large-scale CIB fluctuations are contributed 
by common X-ray sources, which make up ~15%-25% of 
CXB fluctuations, and 

2. ~15%-25% of the large scale CIB fluctuations are con- 
tributed by common X-ray sources, which make up 100% 
of CXB fluctuations. 

Note that y/C(q) is a scale-dependent value. It is important 
to stress that the term “common sources” does not necessarily 
imply that the corresponding parts of the CIB and CXB are 
produced by the same physical sources emitting at both IR and 
X-rays. The IR and X-ray emitters may simply be separated 
by an angle smaller than the Chandra Gaussian beam of ~10" 
corresponding to a physical scale of ~0.1 lr 1 Mpc at z = 1. 
This defines the scale of the individual “objects” in the analysis 
that follows. 

C13 also note that the unresolved CXB fluctuations may be 
contaminated by ionized gas in the Milky Way. Unfortunately, 
this component is difficult to model and subtract but it is not 
expected to exhibit positive cross-correlation with the Galactic 
cirrus and should rather anticorrelate with infrared emitting dust 
clouds. Therefore, the observed X-ray fluctuations are an upper 
limit for the extragalactic CXB component, and its coherence 
with the CIB, quoted atC ~ (0. 15-0. 25 ) 2 , should be considered 
a lower limit. 

4. SOURCES OF THE CIB 

The CIB levels from undetected populations implied by 
the source-subtracted fluctuations require >0.5 nW m -2 sr 1 
(Kashlinsky et al. 2007b) on top of the extrapolated flux from 
known galaxies (Keenan et al. 2010; Ashby et al. 2013). This 
level of CIB is therefore easily accommodated by both direct 
and indirect measurements but the sources of this component 
have not been conclusively identified. However, valuable insight 
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can be obtained from population studies and deep observations 
at other wavelengths. 

4.1. Galaxies 

The extrapolation of faint galaxy populations suggests an 
unresolved CIB of ~0. 1-0.3 nW m~ 2 sr~* at 3.6 /im which is 
mostly produced in the 1 < z < 4 range (Helgason et al. 2012). 
The large-scale clustering of these populations, SF/F < 0.05, 
was found to be insufficient to account for the observed 
CIB fluctuations at large scales (see Figure 6). However, 
provided that unresolved galaxies dominate the unresolved CXB 
fluctuations, they could, in principle, produce enough coherence 
with the CXB to reach C ~ 0.02-0.05 while remaining an 
underdominant component in the CIB fluctuations. In this paper, 
we use the empirical calculation of Helgason et al. (2012) 
in which the CIB from galaxy populations was reconstructed 
using over 230 observed multiwavelength luminosity functions 
(LFs) measured at z < 5. The faint-end regime of this original 
reconstruction was extrapolated to low luminosities, and its 
accuracy has since been verified in deep near-IR source counts 
from the SEDS survey reaching 26 AB mag Ashby et al. (2013). 
We address the X-ray emission from these unresolved galaxies 
in Section 5.1. 

4.2. AGNs 

The AGNs fraction of the resolved CIB sources is small, 
<8%-10% (Treister et al. 2006), but their leading role for 
the CXB makes them important for the interpretation of the 
CIB xCXB correlation. Here, we estimate the CIB production 
from AGNs by constructing an AGN population model from 
7-band LFs of Assef et al. (2011). This sample consists of 
1838 AGNs (Type 1 and 2) at 0 < z < 6 selected in both 
IR and X-ray and is therefore less affected by incompleteness 
and biases seen in purely optical or infrared selected samples. 
The choice of rest-frame 7-band (1.25 /xm) LFs minimizes the 
uncertainty in the k correction, as distant populations observed at 
3. 6-4. 5 /xm emit light at rest frames 0.7 < X <4.5/xmforz < 4. 
We use the pure luminosity evolution parameterization of the 
AGN LF given by Assef et al. (201 1), which includes both host 
galaxy emission and reddening, extrapolating the parameterized 
evolution beyond z > 5 (assuming LDDE model instead does not 
affect our results). The inclusion of the host galaxy is important 
since the X-ray faint AGN tend to be optically obscured and 
dominated by their host galaxy light, which also contributes to 
the CIB. At every distance, the 1.25 /xm emission is corrected 
to l em = 7 0 bs/(l+z) using the average low-resolution AGN 
spectral template of Assef et al. (2010), where 7 0 b s is either 
3.6 and 4.5 /xm for this paper (see also Richards et al. 2006). 
Figure 1 shows a reasonable agreement between this projected 
population and Spitzer/IRAC AGN counts of Treister et al. 
(2006). Inaccuracies in the shape of the counts may arise from 
the fact that our mean AGN spectrum is corrected for host galaxy 
contribution and reddening. 

The typical faint-end slope measured for AGN LFs ranges 
from — 1 .6 to — 1 .3, showing marginal evidence for flattening at 
higher redshifts (Croom et al. 2009; Glikman et al. 20 1 1 ). In our 
model, we conservatively extrapolate a nonevolving faint-end 
slope of p — —1.4 to account for unresolved AGNs but this 
likely overestimates the AGN contribution for several reasons. 
First, we use LFs uncorrected for host galaxy contribution 
that bias the faint-end slope and increase the number of faint 
objects (Hopkins et al. 2007). This effect can be particularly 


S (/xJy) 


1 10 “' 10~ 2 10“ 3 10“ 4 



Figure 1 . Cumulative source number counts of galaxies (green line) compared 
with AGNs (blue line). The models for AGNs are constructed from /-band 
LFs and are uncorrected for host galaxy contribution (explained in the text). 
Corresponding AGN counts data is from Treister et al. (2006, circles). The 
dotted lines show the cases in which the faint-end slope of the AGN LF is 
extrapolated with a constant slope of — 1 .6 (upper) and —1.2 (lower). The green 
line shows the galaxy reconstruction of Helgason et al. (2012) compared with 
SEDS data (Ashby et al. 2013, diamonds). 

(A color version of this figure is available in the online journal.) 

pronounced in the near-IR, as the ratio of host to AGN in 
unobscured objects has a typical maximum at 1.6 /xm. Second, 
the 0.5-2 keV versus 3.6 /xm flux ratio turns toward X/O < 0 
at faint fluxes, suggesting a decreasing importance of AGN 
contribution (Civano et al. 2012). This is consistent with a 
significant number of optically normal galaxies that are seen as 
hosts of low-luminosity X-ray AGNs, (e.g., Barger et al. 2005). 
Diminishing nuclear activity makes the host galaxy dominate at 
low luminosities, and the distinction between galaxy and AGN 
becomes less meaningful. Then, these sources should be largely 
accounted for by our treatment of IR galaxies (see previous 
subsection). In addition, our treatment of normal galaxies (see 
Section 4.1) is based on a compilation of observed galaxy LFs 
which do not typically exclude AGNs. 

From Figure 1, it is clear that normal galaxies are far 
more numerous than AGNs at all flux limits regardless of the 
extrapolation to faint levels. The CIB contribution of AGNs 
is always <10% of the galaxy contribution, and their sky 
surface density is ~2-3 orders of magnitude smaller than the 
~1 arcsec” 2 required to explain the measured CIB fluctuations 
(Kashlinsky et al. 2007b). AGNs would need to live in low- 
mass halos to exist in sufficient numbers at z < 4, and our mask 
eliminates most halos >10 12 M Q where AGNs are typically 
found. Furthermore, a recognizable signature of faint AGNs is 
that their near-IR spectrum should increase with wavelength 
due to their dusty torus emission but this is inconsistent with the 
blue colors of the measured CIB fluctuations (Matsumoto et al. 
2011). Therefore, we assert that the IR emission from AGNs 
themselves is insufficient to produce significant unresolved CIB 
fluctuations. However, this does not necessarily eliminate AGNs 
as sources of the CIB x CXB correlation, as their X-ray emission 
can produce stronger correlation with other IR populations 
sharing common large-scale structures (see Section 7). 

4.3. Diffuse Emission, Intrahalo Light 

Recently, some modeling of the origin of the CIB fluctuations 
has focused on a form of “missing light” associated with galaxy 
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populations but distributed in diffuse structures around masked 
sources. There are several empirical lines of observational 
evidence that argue strongly against such an origin. 

1. There are no spatial correlations between the source- 
subtracted Spitzer / IRAC maps and galaxies detected 
in deep Hubble / ACS maps (0.5-0. 9 /im) < 28 mag 
(Kashlinsky et al. 2007c). However, there are very sig- 
nificant correlations between the ACS galaxies and the un- 
masked Spitzer maps. This means that ACS galaxies and 
any associated diffuse emission, cannot contribute signifi- 
cantly to the large-scale CIB fluctuations found in source- 
subtracted Spitzer data. 

2. Arendt et al. (2010) show that the large-scale CIB fluctua- 
tions are not sensitive to increasing/decreasing the size of 
masked regions around resolved galaxies. Indeed, Figure 17 
of Arendt et al. (2010) shows that there is little variation 
in the CIB fluctuation as the source masking is eroded or 
dilated to masking fractions varying from ~7% to ~46%. 

3. Kashlinsky et al. (2005) and Arendt et al. (2010) show that 
there are no correlations between the source-subtracted CIB 
fluctuations and the identified removed extended sources. 
Moreover, Arendt et al. (2010) constructed artificial halos 
around masked sources and demonstrated that the diffuse 
emission in the final image does not correlate spatially with 
the halos around masked sources that mimic missing light. 

These arguments contradict scenarios invoking any form of 
“missing light” associated with masked galaxies. 

One such scenario was proposed by Cooray et al. (2012b), 
who considered diffuse starlight scattered around and between 
galaxies at z ~ 1-4 as an alternative explanation for the origin 
of the unresolved CIB fluctuations. 8 This can arise from stars 
stripped in mergers or ejected via other processes. In this paper, 
we consider such a diffuse component following definitions in 
Cooray et al. (2012b), deriving the CIB production history as 


d CjHL 
dz 



dn dt , 

LmuAM'Z)— dM — (l+z)" 1 , 
dM dz 


( 2 ) 


where X' = k 0 b s (l +z) using a spectral template of a typical ellip- 
tical galaxy containing old stellar populations from Starburst99 
(Leitherer et al. 1999). This template assumes a 900 Myr old 
stellar population forming at the same time with a Salpeter initial 
mass function (IMF) (1-100 M 0 ) and metallicity Z = 0.008. 
In Figure 4, we show that different choices of spectral energy 
distribution (SED) parameters do not affect the intrahalo light 
(IHL) flux to a great extent (unless the population is very young 
<20 Myr). In all other respects, we follow the formulation 
in Cooray et al. (2012b) recovering a flux of 1 nW m -2 sr -1 
at 3.6 pm which is consistent with but slightly higher than 
the 0.75 nW m -2 sr -1 quoted in Cooray et al. (2012b). We 
are unable to get the large-scale clustering up to the quoted 
SF/F = 10%-15% and, as a result, our large-scale fluctuations 
are a factor of ~2-3 lower; this actually provides a somewhat 
better fit to the measured signal 4.5 pm (see Figure 6). 

As is shown in Figure 4, the CIB from this component. 
Equation (2), exceeds the total CIB from all galaxies already 
at z > 2 according to the reconstruction of Helgason et al. 
(2012). Although this by itself makes the model nonviable, we 
calculate the coherence levels between CIB and diffuse X-ray 
emission for this component. 


8 The empirical results 1-3 above were omitted from the discussion by 
Cooray et al. (2012b). 


5. SOURCES OF THE CXB 
5.1. Normal X-Ray Galaxies 

The bulk of the CXB (~80%-90%) has been resolved into 
point sources, most of which is contributed by AGNs (Hickox 
& Markevitch 2006; Lehmer et al. 2012). However, the number 
counts from the Chandra deep fields reveal that the contribution 
of normal galaxies (mostly XRBs) approaches that of AGNs 
at the faintest levels (Xue et al. 2011; Lehmer et al. 2012). 
A simple extrapolation of the slope of galaxies implies that 
they will ultimately dominate AGNs at <10' 17 erg s -1 cm' 2 
(see Figure 5); this is in fact suggested by deep stacking 
analyses (Cowie et al. 2012). Although Helgason et al. (2012) 
demonstrated the low contribution of known galaxy populations 
to the unresolved CIB, their CXB x CIB amplitude ultimately 
depends on their X-ray properties. 

The X-ray galaxy luminosity function (XLF) derived from 
Chandra and XMM-Newton data is limited to small samples 
of local galaxies of L x > 10 40 erg s' 1 (Norman et al. 2004; 
Tzanavaris & Georgantopoulos 2008) with luminosity evolution 
consistent with L’ oc (1 + z) 23 out to z ~ 1. The X-ray 
emission in galaxies is dominated by a population of compact 
objects accreting from a stellar companion, although hot gas 
can contribute substantially to the soft X-ray flux (<1 keV). To 
gain a better understanding of the X-ray galaxy population and 
its evolution, a popular approach takes advantage of empirical 
correlations of X-ray luminosity with various galaxy properties 
derived from longer wavelength data. The XLF can be related 
to LFs measured at other wavelengths following Avni & 
Tananbaum (1986) 

/ OO 

0(logL Y )PilogL x llogL Y )dlogL Y (3) 

-OO 

(in num Mpc~ 3 dex -1 ), where Y represents a rest-frame band 
which correlates with X-ray luminosity via some specified 
L x -L y relation, and Fling L x \ log Ly) describes the probability 
of a source of L Y having an X-ray luminosity of L x . Here, we 
assume the probability distribution to be Gaussian: 


1 

B(log L x | log L y ) = — exp 

\Z2tti7 


(log L x (L y ) — log L x ) 2 
2a 2 


(4) 

where L X (L Y ) is the X-ray luminosity predicted by a L X -L Y 
relation, and a is the standard deviation of the scatter in the 
measured correlation. 

We consider the X-ray luminosity of normal galaxies to be the 
sum of contributions from high-mass X-ray binaries (HMXBs), 
low-mass X-ray binaries (LMXBs), and hot gas 


Lf = L x (HMXB) + Lx(LMXB) + L z (gas). (5) 


Other forms of galaxy-wide X-ray emission, such as from white 
dwarf binaries and supernova remnants, have been found to be at 
least 10 times weaker, so we ignore their contribution in this pa- 
per (Boroson et al. 2011). A strong correlation is found between 
the SFR and the overall X-ray emission of the galaxy which 
is attributed to active star forming regions producing bright, 
short-lived (<100 Myr) HMXBs. On the other hand, long-lived 
(> 1 Gyr) LMXBs have been found to correlate well with the 
net stellar mass in galaxies. For a full census of X-ray galaxies 
and their different emission mechanism, it is helpful to decom- 
pose the population into late (active) and early (quiescent)-type 
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galaxies, ® tot = ® late + ® early . In the picture that follows, late- 
types can be assumed to contain HMXBs, LMXBs, and hot gas, 
whereas early-types only contain LMXBs and hot gas. 

5.1.1. Late-types 

The tight correlation between X-ray luminosity and SFR 
has established HMXBs as the dominant mechanism for X-ray 
production in late-type galaxies (e.g., Ranalli et al. 2003). For 
the purpose of reconstructing the evolving population of X-ray 
galaxies, we rely on the total infrared LF (8-1000 /im). 9 This 
choice is motivated by several points: (1) reradiated dust 
emission in the FIR is an established tracer of star formation that 
bypasses uncertainties in UV tracers due to obscuration, (2) the 
Lx-Lfir correlation is both very significant and well calibrated, 
and (3) the evolution of the FIR LF has been probed out to z ~ 3 
allowing the estimation of HMXB activity for four-fifths of the 
cosmic time. Since the bulk CXB from galaxies comes from z 
< 3, it is not necessary to rely on UV LFs probing z > 3. In this 
paper, we adopt the FIR LF measured by Magnelli et al. (2009, 
2011) in the 0 < z < 2.5 range from deep Spitzer GOODS/ 
FIDEL data. The measurement is in good agreement with other 
FIR LFs in the literature (Takeuchi et al. 2003; Rodighiero et al. 
2010 ) and is described by a double power law parameterized by 
the characteristic luminosity L*, normalization </>*, and bright 
and faint-end slopes a and /L The evolution is consistent with 
pure luminosity evolution L* oc (1 + z) 36 put to z = 1 with 
very mild evolution of both L * and <fi* at z > 1. Beyond z = 3, 
we assume L* gets fainter with z at the rate implied by the UV 
LF of Bouwens et al. (2011) but this has little impact on our 
results. 

Lehmer et al. (2010) study the Lx-SFR relation for a 
population of star-forming galaxies covering roughly four orders 
of magnitude in SFR, —1.5 < log SFR < 2.5 in M 0 yr -1 . 
The relation is derived at 2-10 keV, where hot gas becomes 
negligible and the emission is predominantly driven by HMXBs 
and LMXBs. They find a local best fit to be 


relation (Equation (6)) slightly overproduces the faintest Chan- 
dra Deep Field-South (CDF-S) counts Lehmer et al. (2012), 
which could also be traced to the evolution of our FIR LF 
which is somewhat steeper than some other measurements in 
the literature (Magnelli et al. 2009). Therefore, we include a 
term of log Lx oc 0.5(1 + 7) in Equation (6), which results 
in good agreement with the data. If we use the observed evolu- 
tion regardless (overproducing faint X-ray galaxies), this gives 
a CIB x CXB signal which is within a factor of 1 .5 larger but our 
final conclusions are unchanged. 

5.1.2. Early-types 

The X-ray emission from quiescent galaxies is thought 
to be dominated by long-lived LMXBs leftover from earlier 
episodes of star formation and hot gas in extended halos. The 
X-ray luminosity of early-types is found to correlate well with 
X-band luminosity, the preferred indicator of stellar mass. For 
the template LF for early-type galaxies, we have chosen the 
local X-band LF of Kochanek et al. (2001) (Two Micron All 
Sky Survey) for z < 0.05 and the evolving LF Arnouts et al. 
(2007) from combined SWIRE-VVDS-CFHTLS data reaching 
Z = 2. These measurements separate the contribution from the 
early-types and late-type to the total LF. The bright (high-mass) 
end of the LF (stellar mass function) is dominated by early- 
types, whereas late-types are much more numerous at the faint 
(low-mass) end. All of the X-band LFs are well described by 
the Schechter function, with parameters L*, </>*, and a measured 
out to z = 2. We fit the evolving parameters in the 0 < z < 2 
range using the functional forms in Helgason et al. (2012) and 
extrapolate them beyond z > 2. 

Boroson et al. (2011) studied the distinct components of 
X-ray emission in nearby early-types in great detail, resolving 
the individual XRBs. The total LMXB luminosity is found to 
correlate with X-band luminosity via the relation 

log L 0 . 3 — 8keV = log L K + 29, (7) 


log L 2 — lOkeV 


0.94 log Lpm + 30. 17, log Lfir ^ 9.6 
0.74 log Lfir + 32.09, log Lfir 9-6 


with a la scatter of roughly 0.4 dex. 10 * It is important to note 
that, because Lehmer et al. (2010) derive the SFR solely based 
on IR luminosity, our model does not depend on an Lfir-SFR 
calibration. The duality of the relationship in Equation (6), over a 
wide range of SFR (0.01-100 M e yr -1 ), arises because LMXBs 
provide a non-negligible contribution to the low SFR regime. 
Indeed, Lehmer et al. (2010) use the X-band luminosity to trace 
stellar mass and find a relationship M* oc SFR 1 1 for SFR < 
5 Mq yr -1 and M* oc SFR 0 3 for SFR >5 M 0 yr -1 , which 
combines SFR and M * in a more physically justified relation, 
Lx = « Af *+/JSFR. Our adopted relation (Equation (6)) therefore 
accounts for both the HMXBs and LMXBs contribution in late- 
type galaxies. We assume an average spectral index of T = 1.8 
to convert to the 0.5-2 keV band. 

Recent evidence seems to indicate that the local Lx-SFR 
relation is not constant with redshift. Basu-Zych et al. 
(2013) find an evolution Lx oc (1 + z) 09 for galaxies with 
SFR > 5 M q yr -1 . However, including this evolution in our 


9 We denote the total infrared (8-1000 fim) as “FIR” to avoid confusing with 
“IR”, which we use to refer to the near-IR. 

10 We have replaced SFR with Lfir (in solar units) based on the original 

relation used by Lehmer et al. (2010), SFR = 9.8 x 10 -11 Lfir. 


where Lk is in Lq,k = 4.82 x 10 32 erg s , and the la scatter is 
» 0.3 dex. We convert the 0.3-8 keV luminosity to the 0.5-2 keV 
band assuming a spectral index of T = 1.8. 

5.1.3. Hot Gas 

It is well known that galaxies contain extended hot halos 
of gas heated above the virial temperature emitting in lines and 
thermal continuum. We find that a luminosity comparable to that 
of XRBs is needed from hot gas in early types to account for the 
number of bright 0.5-2 keV sources (Figure 2), whereas XRBs 
are sufficient to explain the entire late-type population. This is 
consistent with the fact that hot gas in early-types is found to 
constitute a much greater fraction of the total Lx than in late- 
types (Anderson et al. 2013). The unresolved CXB, however, 
cannot contain significant contributions from hot gas for several 
reasons. First, groups and clusters closeby with kT > 1 keV 
are easily detected and removed in the Chandra maps, and the 
joint IR/X-ray mask used in C13 further eliminates galaxies 
residing in >10 12 M 0 (~0.1 keV) halos out to z ~ 2, where 
high-mass systems become increasingly rare. Second, normal 
galaxies have characteristic temperatures of < 1 keV, where their 
spectrum decreases exponentially. At z > 1, their contribution 
quickly redshifts out of the 0.5-2 keV band. Both of these 
considerations act to reduce hot gas contribution in the faint 
unresolved regime, which is dominated by low-mass systems at 
increasingly high redshifts. 
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Figure 2. Cumulative X-ray source counts in soft (left) and hard (right) X-ray bands. The panels show observed galaxy counts from the CDF-S 4 Ms survey (Lehmer 
et al. 2012). The black circles represent all galaxies and the subcontributions from late- and early-type galaxies are shown as purple diamonds and orange triangles, 
respectively. The solid lines show the prediction of our population model in the same color scheme (black, purple, and orange for total, late-, and early-types, 
respectively). We show the extrapolation all the way down to 10 1 erg s -1 cm -2 to emphasize the expected behavior in the unresolved CXB regime. 

(A color version of this figure is available in the online journal.) 


To demonstrate this, we use the hot gas properties predicted 
in semi-analytic models (SAMs) of Guo et al. (2011) mapped 
onto the Millennium Simulations. The hot gas mass is calculated 
based on the baryon content, cooling, and infall rate onto the halo 
(see Guo et al. 2011 for details). We assume a density profile of 
oc (1 + (r/r 0 ) 2 ) _3/3/2 with a constant r 0 = R vir /10 and f> = 2/3 
(Anderson et al. 2013) and apply k corrections assuming a 
thermal continuum spectrum of oc T vl[ exp (— E / kT vn ). In this 
description, the observed [0.5-2 keV] counts are reproduced 
when all systems with masses below 10 13 M 0 are included, 
which is roughly the mass scale of the most massive galaxies. 
However, the hot gas counts flatten toward lower fluxes where 
they are dominated by the HMXB component. This is shown as 
a red dotted line in Figure 5. In Section 5.3 we consider diffuse 
X-ray emission in more complex environments such as groups 
and filaments, including the warm-hot intergalactic medium 
(WHIM). 


5.1.4. The Unresolved CXB from Galaxies 
The cumulative source density seen on the sky is 

r r°° d 3 N 

m>S> = JL dSdlda dSdz ’ 

where the differential is expressed in terms of the XLF as 

d^N 4t rdj(z) dV 

= — 0(Ly,z , 

dSdzdD. K(z) dz 


( 8 ) 


(9) 


where dV /dz = cd\{ 1 + z)~ l dt /dz is the comoving volume 
element and K (z.) is the k-correction. For a power-law spectrum, 
E r , with a spectral index T. the k-correction becomes (1+z) 2 *. 
The power-law slope of X-ray galaxy spectra has been found 
to lie in the 1 < T < 3 range with a mean of ~ 1 .8, which is 
what we assume for the HMXB and LMXB contribution. For 
the hot gas component, we k-correct using a thermal continuum 
spectrum, a T -0,5 exp (—E/kT). 


Figure 2 shows the predicted cumulative X-ray counts of our 
galaxy population models (Equation (8)) and compares them 
with the observed deep counts of Lehmer et al. (2012) showing 
both the early- and late-type contribution. The bright counts rise 
with a near-Euclidian slope and start turning over toward the 
faint regime as expected from both the LF turnover and cosmic 
expansion. At the flux limit of today’s deepest measurements, 
~10 -17 erg s 1 cirT 2 , the source counts are dominated by kp |R 
galaxies (at the knee of the LF) gradually turning over to the 
faint-end regime, which only becomes relevant at much fainter 
levels (<10~ 19 erg s -1 cm -2 ). As long as the flux converges, 
the unresolved CXB fluctuations are not very sensitive to the 
faint-end slope of the XLF. 

The spatial CIBxCXB correlation signal is detected by 
C13 after masking sources resolved in both the Chandra 
and Spitzer/IRAC maps. The common IR/X-ray mask leaves 
~68% of the pixels in the overlapping 8' x 45' region for the 
Fourier analysis. The superior resolution and depth reached in 
the Spitzer exposure maps causes faint X-ray sources to be 
subtracted well below the detection threshold of the Chandra 
maps. The depth of the joint CXBxCIB analysis is therefore 
mostly determined by the IR source subtraction. In this paper, 
we assume a fixed magnitude limit in the Spitzer /\RAC maps 
with X-ray sources removed according to a selection function 
ri{Sx\in\\m), where m\ im refers to the IRAC magnitude limit of 
the near-IR mask. We also tried relaxing the assumption of 
a fixed ni\[ m and instead used the complement of the source 
selection completeness in the SEDS EGS field. This has little 
effect on our results. In order to obtain r]{Sx\m\m/) at flux 
levels inaccessible to current X-ray observatories, we look at 
the distribution of fx //ir predicted in the SAM of Guo et al. 
(2011) mapped onto the Millennium Simulation (for a detailed 
description, we refer to Guo et al. 2011). Mock lightcones based 
on this model were constructed by Henriques et al. (2012). 
Although mostly consistent with galaxy counts in the optical, 
the SAM tends to overestimate the number of small systems 
causing the abundance of faint galaxies to overpredict observed 
3.6 pm and 4.5 pm counts. We therefore apply a post-correction 
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Figure 3. Top: the 4.5 pm vs. 0.5-2 keV flux distribution of normal galaxies 
from the Millennium Simulations lightcones of Henriques et al. (2012). The 
flux limit of the Spitzer / IRAC maps of Kashlinsky et al. (2012) is shown 
as the horizontal line. The color scheme corresponds to the log-scaled flux 
contribution to the total (resolved+unresolved) CIBxCXB SdN/dS, where 
red/blue represents a large/low contribution. Bottom: the unresolved selection 
for X-ray galaxies using an IR threshold of 25 mag extracted from this catalog 
(see Equation (11)). 

(A color version of this figure is available in the online journal.) 

to the galaxy population by shifting excess sources of magnitude 
m by a factor 



Figure 4. Top: the unresolved CIB (4.5 pm) production rate from unresolved 
galaxies, all galaxies, and intrahalo light (green, black, and red). The light 
shaded region shows the range of IHL flux for SED templates with different 
ages (20-900 Myr), metallicities (0.04-0.001 Z Q ), and IMF. Model details are 
explained in Section 4. The AGN contribution falls slightly below the plot range. 
Bottom: the unresolved 0.5-2 keV CXB production histories from galaxies, 
AGNs, and diffuse gas (green, blue, and red). Models are explained in Section 5. 
(A color version of this figure is available in the online journal.) 


. \m(n obs ) - m(n = n obs )\ 

Am = , 10) 

5 

where m(n) are the Millennium SAM derived counts, and 
« 0 bs are the observed counts. We find that for a modification 
factor of 5 = 1.05, the population of Henriques et al. (2012) 
is brought into a good agreement with the observed counts 
while conserving the total number of systems and their redshift 
distribution. Figure 3 shows the 4.5 /im versus 0.5-2 keV 
brightness distribution of normal galaxies according to the 
Millennium SAM where we have used the SFR and M* to 
predict the X-ray luminosity via L x = oiM * + /CSFR of Lehmer 
et al. (2010) and include a scatter of a = 0.4. This is essentially 
the same relation as the Lx-Lyg in Equation (6) with the 
conversion SFR/Lir = 9.8 x 10 11 . The approximate detection 
limit of the Spitzer maps is shown in Figure 3 as a horizontal 
line and the color scheme scales with flux per solid angle, i.e., 
depicting the contribution to the total CXB x CIB background. 
From Figure 3, it is clear that most of the background light 
is resolved and eliminated in the masking process with a 
diminishing contribution from the remaining unresolved sources 
toward the bottom left. We define the selection function of 
0.5-2 keV source removal as the unresolved galaxy fraction 


, c , , N(S x \m>m hm ) 

b(SxHim) = ... c , (11) 

N(S x \m) 

and display it in the bottom panel in Figure 3. This shows that, 
for an IR limit of 25 mag, 90% of sources are removed at 
~3 x 10~ 18 erg s -1 cm~ 2 which is considerably fainter than the 
flux limit of CDF-S (Xue et al. 2011). Systems identified as 


subhalos in the Millennium catalog were removed together with 
its parent halo, provided it is brighter than m\ un . The subhalos 
have little effect on our results. 

Our XFF model allows us to construct the flux production 
rate per solid angle from undetected galaxies as 


dF x 

dz 


f 


n (S)S 


d 3 N 

dSdzdLl 


dS, 


( 12 ) 


where d' N /dS/dz/dQ. comes from Equation (9)b and F x is 
the X-ray flux per solid angle. Figure 4 shows the history of 
the emitted X-rays from galaxies remaining after removing IR 
sources brighter than 25 mag. Very little CXB remains at z < 0.5 
after source subtraction but rises thereafter and peaks at z ~ 1 , 
close to the peak of the star formation history. 


5.2. X-Ray AGN 

The resolved CXB is dominated by AGN populations, which 
have been studied in detail by Chandra and XMM-Newton out 
to z~ 3. For X-ray AGNs, we use the population model of 
Gilli et al. (2007) based on XFFs and evolution of AGNs. The 
models consider the observed XFFs, k corrections, absorption 
distribution, and spectral shapes of AGNs and return the 
observed X-ray flux distribution at any redshift. The models 
have been shown to adequately reproduce source counts, redshift 
distribution, and intrinsic column densities. Our adopted AGN 
population contains sources in the 0 < z < 8 range with 
a wide range in luminosity, 38 < log(L x /ergs ') < 47, to 
allow for very faint unresolved sources. Column densities are 
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Figure 5. Upper: the selection of the unresolved regime for galaxies (darkest 
gray shade) and AGNs (gray shade). Note that, for a constant miR limit, a 
larger fraction of AGN remain unresolved due to their higher X/O ratio and 
greater dispersion. The complement of the X-ray source selection in AEGIS- 
XD (Goulding et al. 2012) is shown for comparison (lightest gray shade). This 
essentially defines the unresolved regime in C13 for X-ray source removal, i.e., 
without additional IR masking. Lower: the source counts of Lehmer et al. (2012) 
compared with our adopted Gilli et al. (2007) model for AGN (blue) and our 
XRB population model (green). The dotted line shows the hot gas contribution 
from virialized halos. 

(A color version of this figure is available in the online journal.) 


20 < logf/Vn/cnr) < 26. The evolution of the AGN XLF in 
Gilli et al. (2007) is modeled with and without an exponential 
decay at z > 2.7 on top of the extrapolated evolution from 
lower redshift parameterization of Hasinger et al. (2005), i.e., 
</>(L,z) = 4>{L,z o)10~°' 43u_Zo - ) with zo = 2.7. We do not include 
the decrease at high-z, which results in a CIBxCXB signal 
within a factor of two of the case with a decline. In Figure 5, 
we show the counts from our adopted AGN model and compare 
with data. 

The extent to which AGN are removed by the joint IR/X-ray 
mask of C 1 3 is estimated based on data from Civano et al. (2012) 
who provide X/1R flux ratios for 1761 sources in the COSMOS 
survey reaching .So.s-zkeV = 1.9 x 1 0“ 16 erg s' 1 cm” 2 . Toward 
faint fluxes AGN tend to become brighter at 3.6 gm deviating 
from the classic X/O = 0. We fit a linear relation to the 3.6 /tm 
versus 0.5-2 keV distribution, mjR = —1.5 log Sx — 1.7, and 
extrapolate to the faint regime with a large Gaussian dispersion 
of a — 1.5 mag. We then apply m iun = 25 to extract the selection 
q(S ) for the X-ray removal (see Equation (11)). The scatter a = 
1 .5 mag is chosen such that the resulting shot noise / J S X N in the 
X-ray power spectrum matches the data. The selection is shown 
in the top panel in Figure 5 demonstrates the extended tail of 
unresolved X-ray AGN caused by the wide dispersion in their IR 
flux (see Civano et al. 2012). In Figure 4 we show the unresolved 
AGN CXB production rate as function of the redshift. The bulk 
of the CXB flux from undetected AGN comes from z ~ 1 . 

5.3. Diffuse Hot Gas and WHIM 

In Section 5.1, we calculated the CXB contribution of hot 
gas heated within galaxies. However, diffuse gas in groups 
and filaments (including the WHIM) also contributes to the 
CXB and was found to dominate the unresolved 0.5-2 keV 


CXB fluctuations of (Cappelluti et al. 2012b). Scaling relations 
indicate that the X-ray masking of C13 removes galaxy clusters 
and groups down to log (M/M 0 ) = (12.5-13.5) (i.e., kT < 
1 .5 keV; Finoguenov et al. 2007). Thus, only the low luminosity 
(low mass) and warm population of galaxy groups and filaments 
contributes to the unresolved CXB. 

Since this class of objects is difficult to model analytically, 
we describe their properties using a set of mock maps from 
Roncarelli et al. (2012), who used a cosmological hydrodynam- 
ical simulation to define the expected X-ray surface brightness 
due to the large-scale structures. The original hydrodynamical 
simulation (see the details in Tornatore et al. 2010) follows the 
evolution of a comoving volume of 37.5 h l Mpc 3 considering 
gravity, hydrodynamics, radiative cooling, and a set of physi- 
cal processes connected with the baryonic component, among 
which a chemical enrichment recipe that allows to follow the 
evolution of seven different metal species of the intergalactic 
medium (IGM). From its outputs, Roncarelli et al. (2012) sim- 
ulated 20 lightcones with a size of ~0.25 deg 2 each, covering 
the redshift interval 0 < z < 1.5. Each pixel of the maps con- 
tains information about the expected observed spectrum in the 
0. 3-2.0 keV band with an energy resolution of 50 eV. The emis- 
sion coming from the IGM was computed assuming an emis- 
sion from an optically thin, collisionally ionized gas (Apec in 
XSPEC) model and considering the abundances of the different 
metal species provided by the simulation. 

These maps/spectra have been convolved with the Chandra 
response in order to reproduce the effective Chandra count 
rates. Since the CXB data of C13 are masked for galaxy 
clusters, we apply a source masking on the simulated maps. We 
have simulated observations with the actual depth of the C13 
field. We have added an artificial isotropic particle and cosmic 
background according to the levels estimated directly from the 
maps of C13. Random Poisson noise was artificially added to 
the image, and we ran a simple wavelet detection with a signal- 
to-noise ratio threshold of 4. We have then excluded all of the 
regions within which the overall encircled signal from sources 
is above 4 a with respect to the background. The unresolved 
CXB production rate averaged over all the realizations is shown 
in Figure 4. 


6. THE ANGULAR AUTO/CROSS-POWER SPECTRUM 
OF MULTIPLE POPULATIONS 


Cosmic background fluctuations can be described by the 
angular power spectrum, Piq). This can be written as the 
sum of the clustering and shot noise of the underlying source 
populations: 

P t ot(<?) = P(q) + PsN, (13) 

where q is the angular wavenumber in rad” 1 . The first term, rep- 
resenting the clustering, can be related to the three-dimensional 
power spectrum of the underlying sources, P(k,z), by projection 
via the Limber equation (Limber 1953): 



Pij(qd c ', z)dz, 

(14) 


where d c (z) is the comoving distance and // (z) = 

Hq^/Q.m( 1 + z) 3 + Oa. The quantities in the square brackets are 
the unresolved flux production rates constructed in the previ- 
ous sections for different source populations, which are denoted 
by the indices i and j running over our three populations: (1) 
galaxies, (2) AGN, and (3) diffuse emission. The summation 
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Figure 6. Auto power spectra of the unresolved CIB fluctuations at 3.6 fim and 4.5 /x m. Data points are from Kashlinsky et al. (2012). The contribution from known 
galaxies and AGNs are show as green and blue line, respectively (the AGN contribution barely visible above the plot range). The hypothetical IHL contribution is 
shown as a red dashed line. We compare this with the original IHL model from Cooray et al. (2012b; dotted lines). All models are convolved with the IRAC beam 
taken from Kashlinsky et al. (2005) and Arendt et al. (2010). 

(A color version of this figure is available in the online journal.) 


results in six terms, three auto power terms ( i = j), and three 
cross terms ( i ^ j) which represent the coupling of different 
populations that live at the same epochs sharing the same en- 
vironments. The shot noise arises from the fluctuation in the 
number of sources within the instrument beam. It only depends 
on the flux distribution of sources and can be expressed as 


Psn = 



(15) 


where 5i; m is the minimum detected source brightness, and i 
denotes the source population as before. Note that the shot noise, 
sometimes called Poisson term, does not have coupling terms 
as it represents a random process and is uncorrelated between 
different populations. Figures 6 and 7 show the modeled auto 
power spectrum of the angular fluctuations in the CIB and the 
CXB respectively, comparing them with current measurements. 

In the description above, “coupling” terms, (i ^ /) refer to the 
correlation of different populations at the same wavelength, not 
cross-power between two wavelengths. The CIB x CXB cross- 
power spectrum can be written 11 



~dF x ~ 


" dF 1R " 

dz 

i 

dz 


P u (qd c \ z)dz, 


(16) 

where the summation results in nine terms representing all 
combinations of cross-correlated X-ray and IR contributions 
from different source populations (note the different summation 
over j compared to the conditional j ^ i in Equation (14)). The 
cross-power of the shot noise term is 


plK.A _ 

r SN — 


n s 


IR C X 


S LK S 


d 2 N 
d S 1R d S x 


dS 1R dS x , 


(17) 


which is added to P lRX (q). However, this expression requires 
additional knowledge of the dS\^/dSx dependence of each 


1 1 Just as the two-point correlation function is related to the auto power 
spectrum, the cross power spectrum is simply the Fourier transform of the 
cross-correlation function, C nm {6 ) = (8 F n (x)8 F m (x +6)). 



Figure 7. Auto power spectra of the unresolved X-ray background fluctuations 
from different populations at the levels of Cl 3. The contribution of normal 
galaxies (mostly XRBs) is shown in green, AGNs in blue, and hot/ warm gas in 
red. We also display the net coupling term of the three components (orange). 
(A color version of this figure is available in the online journal.) 


population in order to be evaluated. As our model construction 
lacks this information, we instead use 


p iR,x 

SN 


/ 


r d P IR d P x 

a /- SN a /~ SN 

dz dz 



dz. 


( 18 ) 


We test the accuracy of this equation using the Millennium 
lightcones, which give P^’ X directly, and find it to be a good 
approximation to Equation (17) (see Figure 9). 


6.1. Halo Model 

Our description of angular fluctuations requires knowl- 
edge of the power spectrum of luminous sources, Pjj(k) 
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(see Equations (14) and (16)). The distribution of sources in- 
side the same collapsed dark matter halos can be related to 
the ACDM matter density field adopting a halo occupation dis- 
tribution (HOD) within a widely used halo model formalism 
(Cooray & Sheth 2002). In this description, the power spectrum 
of clustering can be approximated as the sum of two terms: 

P u (k) = P^ik) + P 2h {k), (19) 

a one-halo term, P lh , describing the correlated fluctuations 
between sources within the same parent halo; and a two- 
halo term, P 2h , arising from the spatial correlation of two 
sources hosted by separate parent halos. As before, the i ^ j 
terms represent coupling terms between different populations, 
whereas for i = j , the expressions reduce to the more familiar 
form of Cooray & Sheth (2002). For source populations i and j , 
these can be written as 



dn (N?N C j + N’fNfiuWM) + N°N s jU 2 (k\M) 

— dM 

dM njiij 


(20) 


Pf-ik) = 

P lin (k)(fl B 


(21) 

where 

B '=J 

f dn AT 

1 z L b(M)dM 

dM tii 

(22) 


B ' = J 

f dn AT 

( -b(M)u(k\M)d M , 

dM tii 

(23) 



z 

Figure 8. Upper: the approximate lowest halo mass removed by the IR mask 
as a function of redshift obtained from the Millennium SAM of Guo et al. 
(2011). The three lines correspond to 80%, 90%, and 95% of systems in the 
Millennium catalog being removed (top to bottom). The dashed line indicates 
the approximate removal threshold of groups/clusters in the X-ray maps (see 
Erfanianfar et al. 2013) Lower: the large scale bias, i e ff (Equation (30)), from 
galaxies, AGNs, and diffuse components. The solid lines show the biasing for 
a masked density field with the mass selection function in the upper panel, 
whereas the dotted lines show the case of no source subtraction. Notice that at 
low redshift, the density field of undetected systems is underbiased. 

(A color version of this figure is available in the online journal.) 


6.2. Halo Occupation, Bias, and Mass Selection 


where i ^ j represents cross terms, and i = j reduces these 
expression to the more familiar form of Cooray & Sheth (2002). 
The individual quantities are defined as follows. 

1. cln/dM is the evolving halo mass function for which we 
use the formalism of Sheth et al. (2001); 

2. z)} (AT shorthand) is the average halo occupation 
of central sources; 

3. ( AT(M , z )) (AT shorthand) is the average halo occupation 
of satellite sources; 

4. n,(z) is the average number density of population i such 
that 


The HOD has been widely investigated for both galaxies 
and AGNs. We adopt standard parameterizations of N C (M) and 
N S {M), which have been found to match observed data. For 
galaxies, we assume the four parameter description of Zheng 
et al. (2005) 12 : 


AT, = _ 
gal -t 


1 + erf 


log M - log M m 

^log M 


(25) 


AT, = - 
gal 2 


1 + erf 


log M — log 2M n 

Ifiog M 


/ M 
V Afsat 


, (26) 


Hi = j ((AT (M , Z)) + (AT (M , Z))) J^dM- (24) 

5. P u "(k. z) is the linear ACDM power spectrum computed 
using the transfer function of Eisenstein & Hu (1998) and 
the adopted cosmological parameters; 

6. u(k\M) is the normalized Fourier transform of the 
Navarro-Frenk- White (NFW) halo profile (Navarro et al. 
1996; Cooray & Sheth 2002); and 

7. b(M,z) is the linear halo bias adopted from the ellipsoidal 
collapse formalism of Sheth et al. (2001). 

It should be kept in mind that all quantities in 
Equations (19)-(22) are evolving with redshift but we have 
omitted this dependence in the expressions to keep the notation 
tidy. 


where M mm is the minimum halo mass that can host a central 
galaxy, and eri og ;M controls the width of the transition of the step 
from zero to one central galaxy. The satellite term has a cutoff 
mass which is twice as large as the one for central galaxies 
and grows as a power law with a slope of a and is normalized 
by M sat . The amplification of the fluctuations through large- 
scale biasing is most sensitive to the choice of M m ; n , which is 
not well known as galaxies in the lowest mass halos are not 
detected. However, at a given redshift, the source subtraction 
removes the massive halos from the top down and limits the 
range of mass scales where the unresolved fluctuation signal 
arises, n < M < M cut (z), where M cut (z) is the lowest mass 
halo removed at z (see Figure 8, upper). We have adopted the 


12 Nf is a shorthand notation for -)) and is in general redshift 

dependent. 
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following parameters of the HOD model motivated by Sloan 
Digital Sky Survey measurements of Zehavi et al. (2011): 
f7 logM = 0.2, M min = 10 10 M Q , M sat = 5 x 10 10 M o , and 
a = 1, where we have deliberately chosen a low cutoff mass, 
M m allowing low-mass halos hosting galaxies well into the 
unresolved regime. 

The HOD of AGNs is less certain due to low number 
statistics but AGNs seem to be preferentially found in halos 
of ~10 12 ' 5 Mq. The HOD has been measured for X-ray selected 
AGNs at z < 1 (Miyaji et al. 2011; Allevato et al. 2012; 
Richardson et al. 2013) and optical quasars out to z ~ 3 (Shen 
et al. 2013; Richardson et al. 2012). Previous studies suggest 
that, compared to optical quasars, X-ray selected AGNs are more 
strongly clustered and reside in more massive host halos but the 
host halo mass range is insufficiently constrained for a definitive 
conclusion (Cappelluti et al. 2012a). Like galaxies, the central 
AGNs are modeled as a softened step function at M mm but the 
satellites are described by a power law with a low-mass rolloff: 


^AGN — -> 


1 + erf 


log M — log M„ 

Glog M 


(27) 


N s — 
iv AGN — 



(28) 


This description has five free parameters: M m i n . the character- 
istic mass scale of the step where the HOD goes from zero 
to a single AGN per halo, with the transition width controlled 
by <J\ogM- The mass at which a halo contains on average one 
satellite AGN is described by M\ ; a is the power-law index con- 
trolling the steepness of the satellite HOD with increasing host 
mass; and M cut is the mass scale below which the satellite HOD 
decays exponentially. We have chosen parameters obtained in a 
numerical study of Chatterjee et al. (2012) that agree with mea- 
sured values when a selection of Lboi > 10 42 erg s 1 is applied. 
We interpolate the redshift evolution of the Lbd > 10 42 erg s 1 
parameters given in Table 2 of Chatterjee et al. (2012). 

The HOD of our diffuse component is somewhat uncertain 
as it does not describe the same population in X-rays and IR, 
i.e., hot/warm gas as opposed to diffuse starlight. We allow the 
diffuse component to trace the NFW halo profile by considering 
a satellite term only, setting the HOD of central sources to zero. 
For diffuse IR, we adopt the parameters from the IHL model 
of Cooray et al. (2012b), discussed in Section 4.3, and assume 
that this component arises in halos in the ~ 10 9 — 10 12 M G range. 
However, gas does not reach sufficient temperatures in such 
small halos but we neglect this by allowing hot gas to live in 
halos anywhere below the mass limit of detected (and removed) 
groups and clusters identified by Erfanianfar et al. (2013) in the 
EGS field (see dashed line in Figure 8). 

When faced with source-subtracted images, the density field 
is modified in the process of masking the brightest sources which 
live in the most massive, and consequently, most biased halos. 
This effect can be accounted for by knowing the mass-dependent 
luminosity distribution, i.e., the conditional LF. Since we do not 
have this information, we use the semi-analytical model of Guo 
et al. (201 1) mapped onto the Millennium Simulations to explore 
the halo mass dependence of source removal in deep IR maps. 
We eliminate all galaxies brighter than ni\\ m = 25 and construct 
a mass selection function defined as the fraction of unresolved 
systems as a function of host halo mass: 


??(M, z) = 


N(M, z I > ffllim) 
N(M, z ) 


(29) 


We display this function in Figure 8 (upper panel), which shows 
the mass scale at which 80%, 90%, and 95% of the systems 
are removed as a function of redshift. We multiply our galaxy 
HOD (Equation (25)) by this function, thereby subtracting the 
massive halos from the density field (including their satellites). 
Since this function is derived from galaxies, we do not apply this 
mass selection to the AGN and diffuse component and instead 
use the X-ray cluster/group detection limits of Erfanianfar et al. 
(2013) as the upper mass limits. The overall effect of this is 
shown in Figure 8. 

The amplification of fluctuations through large scale biasing 
of the sources follows from Equation (22) in the limit where 
u(k\M) —> 1, or equivalently 


bf(z) = 


-I 


dn Nf + N? 
dM i 


- b(M , z)dM, 


(30) 


where the mass-dependent bias comes from the prescription 
of Sheth et al. (2001). The quantity is shown in Figure 8 
for galaxies, AGNs, and diffuse emission. We also show 
the bias without halo subtraction due to i.e., in the 

absence of source masking. Note how the density field becomes 
underbiased at low redshifts where the mask is most effective. 


6.3. Comparing Fluctuation Models with N-body Simulations 

In real measurements, the angular power spectrum is obtained 
directly from the masked and Fourier transformed image, 
P mn (q) = (A m (q)A*(q)) . However, our calculation of P{q) 
relies on the projection of emitting populations via the Limber 
equation with empirically motivated assumptions for their 
HOD. A more sophisticated treatment would link the source 
luminosities to host halo masses in a conditional LF, 0(L(/W)) 
(see, e.g., Cooray 2006; Bethermin et al. 2013). In order to 
test the validity of the approximations made, we make use 
of the Millennium Simulation SAM of Guo et al. (2011) to 
derive the unresolved CIB x CXB power spectrum from galaxy 
populations. The evolving simulation box has been projected to 
construct lightcones which provide 2 deg 2 mock images of the 
extragalactic sky based on the SAM (Henriques et al. 2012). 
From the mock catalogs, we remove all galaxies brighter than 
IRAC 3 . 6 . 4.5 magnitude 25 AB, including substructure associated 
with the parent halo. The X-ray emission is calculated using the 
relation of Lehmer et al. (2010): L x = aM* + ^SFR + (l +z) 0 - 5 , 
where we have added the last term to account for evolution (see 
Section 5.1). 

We calculate the source-subtracted fluctuations directly from 
the mock images, (A m (q)A*(q)) , and compare the results with 
our fluctuation model in Figure 9. Despite the difference in 
approach, this SAM and our population model predict consistent 
fluxes and CIB x CXB cross-power. Caveats worth mentioning 
include the halo resolution limit of Millennium ~10 1() Mq and 
the possibility of spurious power arising from the replication 
of the simulation box required for the lightcone construction 
(see Blaizot et al. 2005). The agreement with our model is 
nevertheless encouraging. 


7. RESULTS 
7.1. CIB Fluctuations 

Figure 6 shows the auto power spectrum of CIB fluctua- 
tions from unresolved galaxies, AGNs, and IHL, comparing 
them with the measurements of C13. The contribution from 
unresolved galaxies (green) is discussed in detail in Helgason 
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Figure 9. Net cross-power spectrum q 2 P lRX (q)/2/7T from normal galaxies in 
units of erg s -1 cm - ’ nW in 2 sr -2 . The green line shows the prediction from 
our population model calculated using the Limber equation and the halo model 
formalism described in the text. The gray areas show the result from directly 
Fourier transforming simulated images, Ar,x(</) = (Air(^)AJ(^)), obtained 
from SAMs based on the Millennium Simulation (Guo et al. 2011). 

(A color version of this figure is available in the online journal.) 

et al. (2012), and the IHL (red) is close to that of Cooray et al. 
(2012b). We only use the default model from Helgason et al. 
(2012) which has been validated in the SEDS survey (Figure 34 
in Ashby et al. 2013), reducing the uncertainties of the faint-end 
extrapolation of the LF. The contribution from AGNs (blue) is 
much smaller due to their low numbers compared to IR galaxies, 
< 10%. In addition, C13 found that their source-subtracted CIB 
power spectrum is independent of the X-ray mask. This means 
that X-ray flux that may be missed by the IR mask, such as from 
the wide wings of the extended Chandra point-spread function 
(PSF), will show up in the CXB power spectrum but will not 
contribute to the cross-power CIB x CXB. 

7.2. CXB Fluctuations 

We find that the CXB power spectrum is dominated by shot 
noise from unresolved AGN (see Figure 7) with the contribution 
from galaxies and gas being considerably lower. However, 
a different study (Cappelluti et al. 2012b) finds the largest 
contribution to come from hot gas. There are various reasons 
for the different findings. First, the X-ray maps of Cappelluti 
et al. (2012b) are deeper (4 Ms) than C13 (800 ks), allowing 
Cappelluti et al. (2012b) to directly mask AGNs to much fainter 
levels. Second, Cappelluti et al. (2012b) modeled AGNs in the 
luminosity range 42 < log (L/L Q ) < 47, whereas we include 
AGNs all the way down to log(L/L 0 ) = 38. This makes a 
substantial difference in the abundance of faint AGNs. 

The net CXB from unresolved galaxies and AGNs is 2.1 x 
10 -13 and 7.9 x 10~ 14 erg s' 1 cirT 2 , respectively. Whereas the 
CXB power spectrum is consistent with being entirely due to 
shot noise from unresolved AGNs, there is a hint of additional 
clustering toward large scales which is not accounted for. No 
reasonable amount of clustering (bias < 30) is sufficient to ac- 
count for enhanced CXB fluctuations on scales >200". If this 
component is real and extragalactic in nature, it could indicate 
the source of the coherence with the CIB fluctuations. How- 
ever, it is important to note that the 0.5-2 keV fluctuations are 
contaminated by foreground emission from the Galaxy, which 
is not sensitive to the removal of extragalactic point sources. 
Therefore, any interpretation of the CXB power spectrum car- 
ries an intrinsic source of uncertainty due to the contribution 
of the Galaxy. A non-negligible Galaxy component at < 1 keV 
could explain why Cl 3 measure a low level of cross-correlation 
between [0.5-2] keV and [2-4.5]-[4.5-7] keV maps. While ir- 
relevant for the CXB x CIB cross-power spectrum, correcting 


for the Galaxy would reduce the measured CXB power spec- 
trum. Additionally, the extended PSF of Chandra could spread 
some fraction of the X-ray point source flux outside the finer IR 
mask. This would not show up in the CIB x CXB cross-power, 
as the large scale CIB fluctuations do not correlate with either 
IR or X-ray removed sources. 

7.3. CIB x CXB Fluctuations 

We start with summarizing the resultant contributions to 
the cross-power from three main components appearing in 
Figure 10: galaxies, AGNs, and diffuse. 

1. Galaxies. The largest contributions to the cross-power 
comes from galaxies (gal IR x galx, green line) and AGNs 
(galiR-AGNx, purple line). A significant galaxy-galaxy 
component is expected because (1) they make up a substan- 
tial unresolved CIB component and (2) dominate AGNs at 
faint X-ray fluxes (see Figure 5). The small-scale cross- 
power is in good agreement with the data, the galiR x galx 
shot noise making up ~30%. However, this fraction de- 
creases toward large scales. In order to explain the shape of 
the CIB x CXB fluctuations at all scales in terms of galax- 
ies only, one needs to alleviate the problem of the low 
clustering with respect to shot noise, as is the case with 
the CIB fluctuations. Simply increasing the flux of the un- 
resolved populations would overproduce the small-scale 
power. A way of increasing the large-scale power with- 
out affecting the shot noise is to enhance the galaxy bias. 
However, this requires bias of >10, which is not expected 
for faint, low-mass systems. We therefore conclude that 
the entire CIB x CXB signal cannot originate from normal 
galaxies unless future measurements show that the cross 
fluctuations (q 2 P(q)/2/it) decrease toward large scales as 
opposed to staying roughly flat, as indicated by the C13 
data. 

2. AGN. Shot noise from unresolved AGNs provides the 
largest contribution, ~60%, to the small-scale CIB x CXB 
power (Figure 10, blue line). This is because a greater frac- 
tion of bright AGNs remains unresolved after IR masking 
(see Figure 5). At large scales, their AGNirxAGNx con- 
tribution is small due to lower flux in both IR and X-rays 
compared with galaxies. However, a substantial contribu- 
tion comes from X-ray AGNs correlating with IR galaxies 
(purple line). As shot noise does not appear in this term, 
the cross-power spectrum has a shape that resembles the 
data but with an amplitude which is more than an order 
of magnitude below the data. To test whether this term 
could be enhanced, we examined the case in which the IR 
source subtraction removes no additional AGNs, i.e., only 
AGNs detected in X-rays are removed. This gives an ampli- 
tude that is one order of magnitude below the data points. 
Enhancing the clustering of the AGN population to the lev- 
els of very biased high-- quasars, corresponding to AGNs 
living in > 1 0 1 3 M 0 halos, still falls below the measured 
levels. In fact, both our AGN population model (Gilli et al. 
2007) and our AGN removal selection (Civano et al. 2012) 
are chosen conservatively and should, if anything, give a 
smaller signal (see Figure 5). 

3. Diffuse. Dispersed starlight around and between masked 
galaxies can share the same environment with diffuse warm 
gas in collapsed halos and filaments. For distant structures, 
however, the thermal spectrum of the ionized gas, ~ 1 keV, 
quickly redshifts out of the 0.5-2 keV band (see Section 5.3) 
and therefore has limited correlations with the IHF that 
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Figure 10. X-ray vs. IR cross-power terms of all source populations compared to measurements from Cappelluti et al. (2013). The total cross-power spectrum 
from all terms is shown as a solid black line. The individual terms are galaxiesjR-galaxiesx (green), galaxiesjR-AGNx (purple), galaxiesiR-diffusex (orange), 
AGNjR-galaxiesx (turquoise), AGNir-AGNx (blue), and diffuseiR-diffusex (red). Auto terms are shown as dashed lines (/ = j), whereas coupling terms (i ■=/=■ j) are 
shown as dashed-dotted lines. 

(A color version of this figure is available in the online journal.) 


mostly arises at different epochs, z ~ 1-4. Despite this, 
the large scale CIB x CXB component arising between IHL 
and warm gas at z < 1 is comparable to that of galaxies and 
AGNs (see Figure 10, red). The diffuse component could 
be made larger if the bulk of IHL emission arose at low-z 
in halos > 1 0 1 2 5 M Q but this would be unlikely to explain 
the entire CIB x CXB data. Coherence between IHL and 
X-ray galaxies /AGNs is also difficult to accommodate for 
the following reasons. In the resolved regime, the point 
sources are masked so any X-ray emission originating in 
the central regions is eliminated with no correlation with 
diffuse IR light outside the mask. In the unresolved regime, 
the IR galaxies themselves should dominate over IHL which 
can only constitute a fraction of the total galaxy light. We 
already account for the coherence of unresolved IR galaxies 
with X-ray galaxies/AGNs. Furthermore, the fact that the 
CIB fluctuations are not sensitive to the X-ray mask argues 
against the removed X-ray sources being responsible for 
the CIB x CXB signal. We note the additional problems 
with the IHL hypothesis in Section 4.3. 

Despite the large uncertainties in the data, there seems to 
be a systematic lack of cross-power at large scales >300", 
where the source clustering is in the linear regime. At small 
scales, our modeled shot noise term is in agreement with 
the data. To better understand these results, we consider a 
simple model composed of linear ACDM clustering and a 
noise term, P(q ) = a\ /Ac dm + « 2 , where a\ and a 2 are free 
parameters and Pacdm is normalized to unity at 1000". This 
model is then convolved with the Chandra response function. 
For 4.5 lira versus [0.5-2 keV], we find best-fit parameters 
a x = (4.3 ± 1.7) x 10- 17 and a 2 = 1.3 ± 0.3 x 1(T 19 in 
erg s -1 cm -2 nW irL 2 sr -1 , resulting in a / 2 / 12 = 1.2 (see 
Figure 1 1). 13 The net power from our model of z < 6 sources 



10 100 1000 
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Figure 11. 4.5 vs. 0.5-2 keV cross-power spectrum. Data points are from Cl 3. 
The dark gray region corresponds to the lcr uncertainty in the best-fit model, 
P(q) = a i Pacdm+^ 2 5 with the two individual terms shown as light gray regions. 
Our model of net contribution from galaxies, AGNs, and diffuse emissions is 
shown as a black dashed line. 

(black line) is a poor fit to the data (x 2 /8 = 2.8), falling 
more than an order of magnitude below the best-fit model. A 
quantitative comparison is shown in Table 2. This distinction 
is not significant at 3.6 /zm versus [0.5-2 keV] due to the 

13 The smallest scale data point at 107 5 is offset with respect to the rest and as 
a result leads to smaller <12 and drives up the x 2 - This is why the gray best-fit 
region seems somewhat below the small-scale data points. If we neglect this 
data point, the best-fit model is brought in perfect agreement with our modeled 
shot noise. 
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Figure 12. Systematic discrepancy of the cross-power spectra of our model and 
the data as a function of angular scale, Pdatat?)/ ^modelG)- The normalized model 
is indicated by the dashed line. The 1 a regions of the best-fit, two-parameter 
model is show in light gray. 


Table 2 

Comparison of our Net Model with the Best-fit Model P(q) = a\ Facdm + a 2 



Clustering 3 

Shot Noise 


(ai x 10 8 * * * * * * * * 17 * ) 

( a 2 x 

10 19 ) 


(3.6 (im) 

(4.5 /zm) 

(3.6 (tm) 

(4.5 /t m) 

Best fit 

2.5 ± 2.1 

4.3 ± 1.7 

1.0 ±0.4 

1.3 ±0.3 

This work 

0.21 

0.11 

1.93 

1.44 


Notes. The power is in units of erg s 1 cm 2 nW m 2 sr 1 . 
a Clustering at 1000". 


large uncertainties but the systematically growing discrepancy 
toward large scales suggests the same behavior as 4.5 /im versus 
[0.5-2 keV]. This is illustrated in Figure 12. 

8. CONCLUSIONS AND DISCUSSION 

In this paper, we have considered known populations of 
X-ray sources at z < 6 in an attempt to explain the measured 

spatial coherence of the unresolved CIB and CXB. The dominant 
contribution to the modeled CIB x CXB signal comes from 
unresolved galaxies containing XRBs and IR galaxies associated 

with X-ray emitting AGNs found within the same large-scale 

structures. However, we find that the combined contribution 

from galaxies, AGNs, and hot gas is unable to produce the large- 

scale cross-power needed to explain the data. At 4.5 fim versus 

0.5-2 keV, the large-scale cross-power is only (2.6% ± 1.0%) 

of the best-fit model but the discrepancy decreases toward small 
angular scales where the shot noise becomes consistent with the 
best fit. At 3.6 /im versus 0.5-2 keV, these identifications are 
not robust due to larger uncertainties in the data. 

Warm gas in <10 13 5 M Q systems can be bright at z < 0.5 
but IR sources are eliminated down to ~10 n Mq in this 
redshift range. A significant source of coherence may arise at 
low z in association with diffuse starlight (IHL). However, the 
correlation of IHL with unresolved X-ray galaxies/AGNs at 
higher redshifts is problematic. Whereas at z = 0, the IHL 
corresponds to ~1% of the total CIB produced by galaxy 
populations, at z > 2, it seems to exceed the total CIB from 
all galaxies (see Figure 4). The lack of correlation between the 

source-subtracted CIB fluctuations and (1) deep Hubble / ACS 
maps (Kashlinsky et al. 2007c) and (2) test halo images (Arendt 
et al. 2010), also makes it difficult to favor missing starlight as 
a dominant component of the fluctuations. 

While other mechanisms capable of producing a correlation 
between X-rays and IR may exist, they are generally expected 


to be much weaker than the dominant forms considered in this 
work: galaxies, AGNs, and hot gas. Thermal emission from 
hot dust ~700 K would have to arise in the local universe as 
it would otherwise redshift out of the near-IR bands and it is 
also inconsistent with the observed blue color of the source- 
subtracted CIB fluctuation in the 2.4 — 4.5 // m range. Because of 
its red colors, any dust dominated component would have to be 
underdominant in the CIB fluctuations while being associated 
with X-ray emission from the dominant CXB component, 
i.e., the Galaxy foreground or obscured AGN. Infrared cirrus 
emission in the Galaxy should absorb X-rays and exhibit a 
negative cross-power contrary to the measurements of Cl 3. In 
the case of obscured AGN, they make up a greater fraction of 
the hard CXB and are less significant in soft X-rays. This is 
not consistent with the fact that C13 detect CIB x CXB cross- 
correlation in the 0.5-2 keV band but not in the harder bands. 
Furthermore, the cross-correlation between the [0.5-2] keV 
band and both [2-4.5] and [4.5-7] keV are small. There is a hint 
of clustering in the CXB power spectrum at >100", which may 
or may not be the source of coherence with the CIB . However, 
the component producing the large scale CIB x CXB cannot 
constitute less than ~15%-20% of the CXB clustering. 

The possibility that the CIB x CXB signal is contributed by 
high-z miniquasars is discussed in Cl 3. Such objects are ex- 
pected to form early and grow rapidly in order to explain the 
population of bright quasars already in place at z ~ 6. Yue et al. 
(2013b) have constructed a population model of highly obscured 
direct collapse black holes that is able to account for (1) net CIB 
measurements and y - ray absorption constraints; (2) the shape 
and amplitude of the source-subtracted CIB fluctuations; (3) the 
unresolved soft CXB level; and (4) the shape and amplitude of 
the spatial coherence in the unresolved CIB x CXB. The inclu- 
sion of such a component improves the CIB x CXB best fit con- 
siderably (see Figure 11). Whether these requirements can be 
realistically satisfied physically by other types of high-z mini- 
quasars, and at the same time stay within limits imposed by 
reionization and black hole mass growth, will be investigated in 
future work. 
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Sciences Fellowship Program Grant — NNX1 1 AO05H. 

REFERENCES 

Ackermann, M., Ajello, M., Allafort, A., et al. 2012, Sci, 338, 1 190 
Allevato, V., Finoguenov, A., Hasinger, G., et al. 2012, ApJ, 758, 47 
Anderson, M. E., Bregman, J. N., & Dai, X. 2013, ApJ, 762, 106 
Arendt, R. G., Kashlinsky, A., Moseley, S. H., & Mather, J. 2010, ApJS, 

186, 10 

Arnouts, S., Walcher, C. J., Le Fevre, O., et al. 2007, A&A, 476, 137 
Ashby, M. L. N., Willner, S. P., Fazio, G. G., et al. 2013, ApJ, 769, 80 
Assef, R. J., Kochanek, C. S., Ashby, M. L. N., et al. 2011, ApJ, 728, 56 
Assef, R. J., Kochanek, C. S., Brodwin, M., et al. 2010, ApJ, 713, 970 
Avni, Y., & Tananbaum, H. 1986, ApJ, 305, 83 
Barger, A. J., Cowie, L. L., Mushotzky, R. F., et al. 2005, AJ, 129, 578 
Basu-Zych, A. R., Lehmer, B. D., Hornschemeier, A. E., et al. 2013, ApJ, 

762, 45 

Bethermin, M., Wang, L., Dore, O., et al. 2013, A&A, 557, A66 
Blaizot, J., Wadadekar, Y., Guiderdoni, B., et al. 2005, MNRAS, 360, 159 
Boroson, B., Kim, D.-W., & Fabbiano, G. 2011, ApJ, 729, 12 
Bouwens, R. J., Illingworth, G. D., Oesch, P. A., et al. 2011, ApJ, 737, 90 
Cappelluti, N., Allevato, V., & Finoguenov, A. 2012a, AdAst, 2012, 853701 
Cappelluti, N., Kashlinsky, A., Arendt, R. G., et al. 2013, ApJ, 769, 68 
Cappelluti, N., Ranalli, P., Roncarelli, M., et al. 2012b, MNRAS, 427, 651 
Chatterjee, S., Degraf, C., Richardson, J., et al. 2012, MNRAS, 419, 2657 


15 


The Astrophysical Journal, 785:38 (16pp), 2014 April 10 


Helgason et al. 


Civano, F., Elvis, M.. Brusa, M., et al. 2012, ApJS, 201, 30 
Cooray, A. 2006, MNRAS, 365, 842 

Cooray, A., Bock, J. J.. Keatin, B., Lange, A. E., & Matsumoto, T. 2004, ApJ, 
606, 611 

Cooray, A., Gong, Y., Smidt, J., & Santos, M. G. 2012a, ApJ, 756, 92 
Cooray, A., & Sheth, R. 2002, PhR, 372, 1 

Cooray, A., Smidt, J., de Bernardis, F., et al. 2012b, Natur, 490, 514 
Cowie, L. L., Barger, A. J., & Hasinger, G. 2012, ApJ, 748, 50 
Croom, S. M„ Richards, G. T„ Shanks, T„ et al. 2009, MNRAS, 399, 1755 
Dwek, E„ & Arendt, R. G. 1998, ApJL, 508, L9 
Eisenstein, D. J., & Hu, W. 1998, ApJ, 496, 605 
Erfanianfar, G., Finoguenov, A., Tanaka, M., et al. 2013, ApJ, 765, 1 17 
Fernandez, E. R„ Komatsu, E., Iliev, I. T., & Shapiro, P. R. 2010, ApJ, 
710, 1089 

Finoguenov, A., Guzzo, L., Hasinger, G., et al. 2007, ApJS, 172, 182 
Gilli, R., Comastri, A., & Hasinger, G. 2007, A&A, 463, 79 
Glikman, E., Djorgovski, S. G., Stern, D., et al. 2011, ApJL, 728, L26 
Goulding, A. D., Forman, W. R., Hickox, R. C., et al. 2012, ApJS, 202, 6 
Guo, Q.Twhite, S., Boylan-Kolchin, M„ et al. 2011, MNRAS, 413, 101 
Hasinger, G., Miyaji, T., & Schmidt, M. 2005, A&A, 441, 417 
Helgason, K., Ricotti, M., & Kashlinsky, A. 2012, ApJ, 752, 1 13 
Henriques, B. M. B., White, S. D. M., Lemson, G., et al. 2012, MNRAS, 421, 
2904 

H.E.S.S. Collaboration. 2013, A&A, 550, A4 

Hickox, R. C., & Markevitch, M. 2006, ApJ, 645, 95 

Hopkins, P. F., Richards, G. T., & Hernquist, L. 2007, ApJ, 654, 731 

Kashlinsky, A. 2005, PhR, 409, 361 

Kashlinsky, A., Arendt, R.. Gardner, J. P., Mather, J. C., & Moseley, S. H. 
2004, ApJ, 608, 1 

Kashlinsky, A., Arendt, R. G., Ashby, M. L. N., et al. 2012, ApJ, 753, 63 
Kashlinsky, A., Arendt, R. G., Mather, J., & Moseley, S. H. 2005, Natur, 
438, 45 

Kashlinsky, A., Arendt, R. G., Mather, J., & Moseley, S. H. 2007a, ApJL, 
654, L5 

Kashlinsky, A., Arendt, R. G., Mather, J., & Moseley, S. H. 2007b, ApJL, 
654, LI 

Kashlinsky, A., Arendt, R. G., Mather, J., & Moseley, S. H. 2007c, ApJL, 
666, LI 

Keenan, R. C., Trouille, L., Barger, A. J., Cowie, L. L., & Wang, W.-H. 
2010, ApJS, 186, 94 

Kochanek, C. S., Pahre, M. A., Falco, E. E., et al. 2001, ApJ, 560, 566 
Lehmer, B. D., Alexander, D. M., Bauer, F. E., et al. 2010, ApJ, 724, 559 


Lehmer, B. D„ Xue, Y. Q„ Brandt, W. N„ et al. 2012, ApJ, 752, 46 
Leitherer, C., Schaerer, D., Goldader, J. D., et al. 1999, ApJS, 123, 3 
Limber, D. N. 1953, ApJ, 1 17, 134 

Magnelli, B„ Elbaz, D„ Chary, R. R„ et al. 2009, A&A, 496, 57 
Magnelli, B„ Elbaz, D„ Chary, R. R„ et al. 2011, A&A, 528, A35 
Matsumoto, T., Matsuura, S., Murakami, H., et al. 2005, ApJ, 626, 31 
Matsumoto, T., Seo, H. J., Jeong, W.-S., et al. 2011, ApJ, 742, 124 
Meyer, M„ Raue, M., Mazin, D., & Homs, D. 2012, A&A, 542, A59 
Miyaji, T., Krumpe, M., Coil, A. L., & Aceves, H. 201 1, ApJ, 726, 83 
Navarro, J. F„ Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 563 
Norman, C., Ptak, A., Homschemeier, A., et al. 2004, ApJ, 607, 721 
Persic, M„ & Rephaeli, Y. 2007, A&A, 463, 481 

Planck Collaboration, Ade, P. A. R., Aghanim, N., Armitage-Caplan, C., et al. 
2013, arXiv: 1303.5076 

Ranalli, P, Comastri, A., & Setti, G. 2003, A&A, 399, 39 
Richards, G. T., Lacy, M., Storrie-Lombardi, L. J., et al. 2006, ApJS, 
166, 470 

Richardson, J., Chatterjee, S., Zheng, Z., Myers, A. D., & Hickox, R. 2013, ApJ, 
774, 143 

Richardson. J., Zheng, Z., Chatterjee, S., Nagai, D., & Shen, Y. 2012, ApJ, 
755, 30 

Rodighiero, G., Vaccari, M., Franceschini, A., et al. 2010, A&A, 5 15, A8 
Roncarelli. M., Cappelluti, N., Borgani, S., Branchini, E., & Moscardini, L. 
2012, MNRAS, 424, 1012 

Salvaterra, R., Magliocchetti, M., Ferrara, A., & Schneider, R. 2006, MNRAS, 
368, L6 

Shen, Y., McBride, C. K., White, M., et al. 2013, ApJ, 778, 98 
Sheth, R. K.. Mo, H. J., & Tormen, G. 2001, MNRAS, 323, 1 
Takeuchi, T. T„ Yoshikawa, K„ & Ishii, T. T. 2003, ApJL, 587, L89 
Thompson, R. I., Eisenstein, D., Fan, X., Rieke, M., & Kennicutt, R. C. 
2007, ApJ, 657, 669 

Tomatore, L., Borgani, S., Viel, M., & Springel, V. 2010, MNRAS, 402, 1911 
Treister, E., Urry, C. M., Van Duyne, J., et al. 2006, ApJ, 640, 603 
Tsumura, K., Matsumoto, T., Matsuura, S., Sakon, I., & Wada, T. 2013, PASJ, 
65, 121 

Tzanavaris, P, & Georgantopoulos, I. 2008, A&A, 480, 663 
Xue, Y. Q., Luo, B„ Brandt, W. N„ et al. 2011, ApJS, 195, 10 
Yue, B., Ferrara, A., Salvaterra, R., & Chen, X. 2013a, MNRAS, 431, 383 
Yue, B., Ferrara, A., Salvaterra, R., Xu, Y., & Chen, X. 2013b, MNRAS, 
433. 1556 

Zehavi, I., Zheng, Z., Weinberg, D. H., et al. 2011, ApJ, 736, 59 
Zheng, Z., Berlind, A. A., Weinberg, D. H., et al. 2005, ApJ, 633, 791 


16 


